Sequential genome-wide association studies for monitoring adverse events in the clinical evaluation of new drugs

Pharmacovigilance, the monitoring of adverse events (AEs), is an integral part in the clinical evaluation of a new drug. Until recently, attempts to relate the incidence of AEs to putative causes have been restricted to the evaluation of simple demographic and environmental factors. The advent of large-scale genotyping, however, provides an opportunity to look for associations between AEs and genetic markers, such as single nucleotides polymorphisms (SNPs). It is envisaged that a very large number of SNPs, possibly over 500 000, will be used in pharmacovigilance in an attempt to identify any genetic difference between patients who have experienced an AE and those who have not. We propose a sequential genome-wide association test for analysing AEs as they arise, allowing evidence-based decision-making at the earliest opportunity. This gives us the capability of quickly establishing whether there is a group of patients at high-risk of an AE based upon their DNA. Our method provides a valid test which takes account of linkage disequilibrium and allows for the sequential nature of the procedure. The method is more powerful than using a correction, such as idák, that assumes that the tests are independent. Copyright © 2006 John Wiley & Sons, Ltd.

A likelihood ratio appropach to family-based association studies with covariates

We introduce a procedure for association based analysis of nuclear families that allows for dichotomous and more general measurements of phenotype and inclusion of covariate information. Standard generalized linear models are used to relate phenotype and its predictors. Our test procedure, based on the likelihood ratio, unifies the estimation of all parameters through the likelihood itself and yields maximum likelihood estimates of the genetic relative risk and interaction parameters. Our method has advantages in modelling the covariate and gene-covariate interaction terms over recently proposed conditional score tests that include covariate information via a two-stage modelling approach. We apply our method in a study of human systemic lupus erythematosus and the C-reactive protein that includes sex as a covariate.

Genetic association analysis of vitamin D pathway with obesity traits

Objective:Observational studies have examined the link between vitamin D deficiency and obesity traits. Some studies have reported associations between vitamin D pathway genes such as VDR, GC and CYP27B1 with body mass index (BMI) and waist circumference (WC); however, the findings have been inconsistent. Therefore, we investigated the involvement of vitamin D metabolic pathway genes in obesity-related traits in a large population-based study.Methods:We undertook a comprehensive analysis between 100 tagging single nucleotide polymorphisms (tagSNPs) in genes encoding for DHCR7, CYP2R1, VDBP, CYP27B1, CYP27A1, CYP24A1, VDR and RXRG, and obesity traits in 5224 participants (aged 45 years) in the 1958 British birth cohort (1958BC). We further extended our analyses to investigate the associations between SNPs and obesity traits using the summary statistics from the GIANT (Genetic Investigation of Anthropometric Traits) consortium (n=123 865).Results:In the 1958BC (n=5224), after Bonferroni correction, none of the tagSNPs were associated with obesity traits except for one tagSNP from CYP24A1 that was associated with waist-hip ratio (WHR) (rs2296239, P=0.001). However, the CYP24A1 SNP was not associated with BMI-adjusted WHR (WHRadj) in the 1958BC (rs2296239, P=1.00) and GIANT results (n=123 865, P=0.18). There was also no evidence for an interaction between the tagSNPs and obesity on BMI, WC, WHR and WHRadj in the 1958BC. In the GIANT consortium, none of the tagSNPs were associated with obesity traits.Conclusions:Despite a very large study, our findings suggest that the vitamin D pathway genes are unlikely to have a major role in obesity-related traits in the general population.

Improving power and accuracy of genome-wide association studies via a multi-locus mixed linear model methodology

Genome-wide association studies (GWAS) have been widely used in genetic dissection of complex traits. However, common methods are all based on a fixed-SNP-effect mixed linear model (MLM) and single marker analysis, such as efficient mixed model analysis (EMMA). These methods require Bonferroni correction for multiple tests, which often is too conservative when the number of markers is extremely large. To address this concern, we proposed a random-SNP-effect MLM (RMLM) and a multi-locus RMLM (MRMLM) for GWAS. The RMLM simply treats the SNP-effect as random, but it allows a modified Bonferroni correction to be used to calculate the threshold p value for significance tests. The MRMLM is a multi-locus model including markers selected from the RMLM method with a less stringent selection criterion. Due to the multi-locus nature, no multiple test correction is needed. Simulation studies show that the MRMLM is more powerful in QTN detection and more accurate in QTN effect estimation than the RMLM, which in turn is more powerful and accurate than the EMMA. To demonstrate the new methods, we analyzed six flowering time related traits in Arabidopsis thaliana and detected more genes than previous reported using the EMMA. Therefore, the MRMLM provides an alternative for multi-locus GWAS.

Candidate genes for obesity-susceptibility show enriched association within a large genome-wide association study for BMI

Before the advent of genome-wide association studies (GWASs), hundreds of candidate genes for obesity-susceptibility had been identified through a variety of approaches. We examined whether those obesity candidate genes are enriched for associations with body mass index (BMI) compared with non-candidate genes by using data from a large-scale GWAS. A thorough literature search identified 547 candidate genes for obesity-susceptibility based on evidence from animal studies, Mendelian syndromes, linkage studies, genetic association studies and expression studies. Genomic regions were defined to include the genes ±10 kb of flanking sequence around candidate and non-candidate genes. We used summary statistics publicly available from the discovery stage of the genome-wide meta-analysis for BMI performed by the genetic investigation of anthropometric traits consortium in 123 564 individuals. Hypergeometric, rank tail-strength and gene-set enrichment analysis tests were used to test for the enrichment of association in candidate compared with non-candidate genes. The hypergeometric test of enrichment was not significant at the 5% P-value quantile (P = 0.35), but was nominally significant at the 25% quantile (P = 0.015). The rank tail-strength and gene-set enrichment tests were nominally significant for the full set of genes and borderline significant for the subset without SNPs at P < 10(-7). Taken together, the observed evidence for enrichment suggests that the candidate gene approach retains some value. However, the degree of enrichment is small despite the extensive number of candidate genes and the large sample size. Studies that focus on candidate genes have only slightly increased chances of detecting associations, and are likely to miss many true effects in non-candidate genes, at least for obesity-related traits.

High throughput, high resolution selection of polymorphic microsatellite loci for multiplex analysis

Background Large-scale genetic profiling, mapping and genetic association studies require access to a series of well-characterised and polymorphic microsatellite markers with distinct and broad allele ranges. Selection of complementary microsatellite markers with non-overlapping allele ranges has historically proved to be a bottleneck in the development of multiplex microsatellite assays. The characterisation process for each microsatellite locus can be laborious and costly given the need for numerous, locus-specific fluorescent primers. Results Here, we describe a simple and inexpensive approach to select useful microsatellite markers. The system is based on the pooling of multiple unlabelled PCR amplicons and their subsequent ligation into a standard cloning vector. A second round of amplification utilising generic labelled primers targeting the vector and unlabelled locus-specific primers targeting the microsatellite flanking region yield allelic profiles that are representative of all individuals contained within the pool. Suitability of various DNA pool sizes was then tested for this purpose. DNA template pools containing between 8 and 96 individuals were assessed for the determination of allele ranges of individual microsatellite markers across a broad population. This helped resolve the balance between using pools that are large enough to allow the detection of many alleles against the risk of including too many individuals in a pool such that rare alleles are over-diluted and so do not appear in the pooled microsatellite profile. Pools of DNA from 12 individuals allowed the reliable detection of all alleles present in the pool. Conclusion The use of generic vector-specific fluorescent primers and unlabelled locus-specific primers provides a high resolution, rapid and inexpensive approach for the selection of highly polymorphic microsatellite loci that possess non-overlapping allele ranges for use in large-scale multiplex assays.

Why are autism spectrum conditions more prevalent in males?

Autism Spectrum Conditions (ASC) are much more common in males, a bias that may offer clues to the etiology of this condition. Although the cause of this bias remains a mystery, we argue that it occurs because ASC is an extreme manifestation of the male brain. The extreme male brain (EMB) theory, first proposed in 1997, is an extension of the Empathizing-Systemizing (E-S) theory of typical sex differences that proposes that females on average have a stronger drive to empathize while males on average have a stronger drive to systemize. In this first major update since 2005, we describe some of the evidence relating to the EMB theory of ASC and consider how typical sex differences in brain structure may be relevant to ASC. One possible biological mechanism to account for the male bias is the effect of fetal testosterone (fT). We also consider alternative biological theories, the X and Y chromosome theories, and the reduced autosomal penetrance theory. None of these theories has yet been fully confirmed or refuted, though the weight of evidence in favor of the fT theory is growing from converging sources (longitudinal amniocentesis studies from pregnancy to age 10 years old, current hormone studies, and genetic association studies of SNPs in the sex steroid pathways). Ultimately, as these theories are not mutually exclusive and ASC is multi-factorial, they may help explain the male prevalence of ASC.

Variation in the autism candidate gene GABRB3 modulates tactile sensitivity in typically developing children

Background: Autism spectrum conditions have a strong genetic component. Atypical sensory sensitivities are one of the core but neglected features of autism spectrum conditions. GABRB3 is a well-characterised candidate gene for autism spectrum conditions. In mice, heterozygous Gabrb3 deletion is associated with increased tactile sensitivity. However, no study has examined if tactile sensitivity is associated with GABRB3 genetic variation in humans. To test this, we conducted two pilot genetic association studies in the general population, analysing two phenotypic measures of tactile sensitivity (a parent-report and a behavioural measure) for association with 43 SNPs in GABRB3. Findings: Across both tactile sensitivity measures, three SNPs (rs11636966, rs8023959 and rs2162241) were nominally associated with both phenotypes, providing a measure of internal validation. Parent-report scores were nominally associated with six SNPs (P <0.05). Behaviourally measured tactile sensitivity was nominally associated with 10 SNPs (three after Bonferroni correction). Conclusions: This is the first human study to show an association between GABRB3 variation and tactile sensitivity. This provides support for the evidence from animal models implicating the role of GABRB3 variation in the atypical sensory sensitivity in autism spectrum conditions. Future research is underway to directly test this association in cases of autism spectrum conditions.

Family-based association analysis with ordered categorical phenotypes, covariates and interactions

Genetic association analyses of family-based studies with ordered categorical phenotypes are often conducted using methods either for quantitative or for binary traits, which can lead to suboptimal analyses. Here we present an alternative likelihood-based method of analysis for single nucleotide polymorphism (SNP) genotypes and ordered categorical phenotypes in nuclear families of any size. Our approach, which extends our previous work for binary phenotypes, permits straightforward inclusion of covariate, gene-gene and gene-covariate interaction terms in the likelihood, incorporates a simple model for ascertainment and allows for family-specific effects in the hypothesis test. Additionally, our method produces interpretable parameter estimates and valid confidence intervals. We assess the proposed method using simulated data, and apply it to a polymorphism in the c-reactive protein (CRP) gene typed in families collected to investigate human systemic lupus erythematosus. By including sex interactions in the analysis, we show that the polymorphism is associated with anti-nuclear autoantibody (ANA) production in females, while there appears to be no effect in males.

Genetic predisposition influences plasma lipids of participants on habitual diet, but not the response to reductions in dietary intake of saturated fatty acids

Objective: SNPs identified from genome wide association studies associate with lipid risk markers of cardiovascular disease. This study investigated whether these SNPs altered the plasma lipid response to diet in the ‘RISCK’ study cohort. Methods: Participants (n = 490) from a dietary intervention to lower saturated fat by replacement with carbohydrate or monounsaturated fat, were genotyped for 39 lipid-associated SNPs. The association of each individual SNP, and of the SNPs combined (using genetic predisposition scores), with plasma lipid concentrations was assessed at baseline, and on change in response to 24 weeks on diets. Results: The associations between SNPs and lipid concentrations were directionally consistent with previous findings. The genetic predisposition scores were associated with higher baseline concentrations of plasma total(P = 0.02) and LDL (P = 0.002) cholesterol, triglycerides (P = 0.001) and apolipoprotein B (P = 0.004), and with lower baseline concentrations of HDL cholesterol (P < 0.001) and apolipoprotein A-I (P < 0.001). None of the SNPs showed significant association with the reduction of plasma lipids in response to the dietary interventions and there was no evidence of diet-gene interactions. Conclusion: Results from this exploratory study have shown that increased genetic predisposition was associated with an unfavourable plasma lipid profile at baseline, but did not influence the improvement in lipid profiles by the low-saturated-fat diets.

Genetic predisposition to type 2 diabetes is associated with impaired insulin secretion but does not modify insulin resistance or secretion in response to an intervention to lower dietary saturated fat

Genome-wide association studies have identified SNPs reproducibly associated with type 2 diabetes (T2D). We examined the effect of genetic predisposition to T2D on insulin sensitivity and secretion using detailed phenotyping in overweight individuals with no diagnosis of T2D. Furthermore, we investigated whether this genetic predisposition modifies the responses in beta-cell function and insulin sensitivity to a 24-week dietary intervention. We genotyped 25 T2D-associated SNPs in 377 white participants from the RISCK study. Participants underwent an IVGTT prior to and following a dietary intervention that aimed to lower saturated fat intake by replacement with monounsaturated fat or carbohydrate. We composed a genetic predisposition score (T2D-GPS) by summing the T2D risk-increasing alleles of the 25 SNPs and tested for association with insulin secretion and sensitivity at baseline, and with the change in response to the dietary intervention. At baseline, a higher T2D-GPS was associated with lower acute insulin secretion (AIRg 4% lower/risk allele, P = 0.006) and lower insulin secretion for a given level of insulin sensitivity, assessed by the disposition index (DI 5% lower/risk allele, P = 0.002), but not with insulin sensitivity (Si). T2D-GPS did not modify changes in insulin secretion, insulin sensitivity or the disposition index in response to the dietary interventions to lower saturated fat. Participants genetically predisposed to T2D have an impaired ability to compensate for peripheral insulin resistance with insulin secretion at baseline, but this does not modify the response to a reduction in dietary saturated fat through iso-energetic replacement with carbohydrate or monounsaturated fat.

Genetic structure and linkage disequilibrium in landrace populations of barley in Sardinia

Multilocus digenic linkage disequilibria (LD) and their population structure were investigated in eleven landrace populations of barley (Hordeum vulgare ssp. vulgare L.) in Sardinia, using 134 dominant simple-sequence amplified polymorphism markers. The analysis of molecular variance for these markers indicated that the populations were partially differentiated (F ST = 0.18), and clustered into three geographic areas. Consistent with this population pattern, STRUCTURE analysis allocated individuals from a bulk of all populations into four genetic groups, and these groups also showed geographic patterns. In agreement with other molecular studies in barley, the general level of LD was low (13 % of locus pairs, with P < 0.01) in the bulk of 337 lines, and decayed steeply with map distance between markers. The partitioning of multilocus associations into various components indicated that genetic drift and founder effects played a major role in determining the overall genetic makeup of the diversity in these landrace populations, but that epistatic homogenising or diversifying selection was also present. Notably, the variance of the disequilibrium component was relatively high, which implies caution in the pooling of barley lines for association studies. Finally, we compared the analyses of multilocus structure in barley landrace populations with parallel analyses in both composite crosses of barley on the one hand and in natural populations of wild barley on the other. Neither of these serves as suitable mimics of landraces in barley, which require their own study. Overall, the results suggest that these populations can be exploited for LD mapping if population structure is controlled.

Genetic variation in GABRB3 is associated with Asperger syndrome and multiple endophenotypes relevant to autism.

BACKGROUND: Autism spectrum conditions (ASC) are associated with deficits in social interaction and communication, alongside repetitive, restricted, and stereotyped behavior. ASC is highly heritable. The gamma-aminobutyric acid (GABA)-ergic system has been associated consistently with atypicalities in autism, in both genetic association and expression studies. A key component of the GABA-ergic system is encoded by the GABRB3 gene, which has been previously implicated both in ASC and in individual differences in empathy. METHODS: In this study, 45 genotyped single nucleotide polymorphisms (SNPs) within GABRB3 were tested for association with Asperger syndrome (AS), and related quantitative traits measured through the following tests: the Empathy Quotient (EQ), the Autism Spectrum Quotient (AQ), the Systemizing Quotient-Revised (SQ-R), the Embedded Figures Test (EFT), the Reading the Mind in the Eyes Test (RMET), and the Mental Rotation Test (MRT). Two-loci, three-loci, four-loci haplotype analyses, and one seven-loci haplotype analysis were also performed in the AS case--control sample. RESULTS: Three SNPs (rs7180158, rs7165604, rs12593579) were significantly associated with AS, and two SNPs (rs9806546, rs11636966) were significantly associated with EQ. Two SNP-SNP pairs, rs12438141-rs1035751 and rs12438141-rs7179514, showed significant association with variation in the EFT scores. One SNP-SNP pair, rs7174437-rs1863455, was significantly associated with variation in the MRT scores. Additionally, a few haplotypes, including a 19 kb genomic region that formed a linkage disequilibrium (LD) block in our sample and contained several nominally significant SNPs, were found to be significantly associated with AS. CONCLUSION: The current study confirms the role of GABRB3 as an important candidate gene in both ASC and normative variation in related endophenotypes.

Genetic variant rs17225178 in the ARNT2 gene is associated with Asperger Syndrome

Background Autism Spectrum Conditions (ASC) are neurodevelopmental conditions characterized by difficulties in communication and social interaction, alongside unusually repetitive behaviours and narrow interests. Asperger Syndrome (AS) is one subgroup of ASC and differs from classic autism in that in AS there is no language or general cognitive delay. Genetic, epigenetic and environmental factors are implicated in ASC and genes involved in neural connectivity and neurodevelopment are good candidates for studying the susceptibility to ASC. The aryl-hydrocarbon receptor nuclear translocator 2 (ARNT2) gene encodes a transcription factor involved in neurodevelopmental processes, neuronal connectivity and cellular responses to hypoxia. A mutation in this gene has been identified in individuals with ASC and single nucleotide polymorphisms (SNPs) have been nominally associated with AS and autistic traits in previous studies. Methods In this study, we tested 34 SNPs in ARNT2 for association with AS in 118 cases and 412 controls of Caucasian origin. P values were adjusted for multiple comparisons, and linkage disequilibrium (LD) among the SNPs analysed was calculated in our sample. Finally, SNP annotation allowed functional and structural analyses of the genetic variants in ARNT2. We tested the replicability of our result using the genome-wide association studies (GWAS) database of the Psychiatric Genomics Consortium (PGC). Results We report statistically significant association of rs17225178 with AS. This SNP modifies transcription factor binding sites and regions that regulate the chromatin state in neural cell lines. It is also included in a LD block in our sample, alongside other genetic variants that alter chromatin regulatory regions in neural cells. Conclusions These findings demonstrate that rs17225178 in the ARNT2 gene is associated with AS and support previous studies that pointed out an involvement of this gene in the predisposition to ASC.

A comprehensive meta-analysis of common genetic variants in autism spectrum conditions

Background Autism spectrum conditions (ASC) are a group of neurodevelopmental conditions characterized by difficulties in social interaction and communication alongside repetitive and stereotyped behaviours. ASC are heritable, and common genetic variants contribute substantial phenotypic variability. More than 600 genes have been implicated in ASC to date. However, a comprehensive investigation of candidate gene association studies in ASC is lacking. Methods In this study, we systematically reviewed the literature for association studies for 552 genes associated with ASC. We identified 58 common genetic variants in 27 genes that have been investigated in three or more independent cohorts and conducted a meta-analysis for 55 of these variants. We investigated publication bias and sensitivity and performed stratified analyses for a subset of these variants. Results We identified 15 variants nominally significant for the mean effect size, 8 of which had P values below a threshold of significance of 0.01. Of these 15 variants, 11 were re-investigated for effect sizes and significance in the larger Psychiatric Genomics Consortium dataset, and none of them were significant. Effect direction for 8 of the 11 variants were concordant between both the datasets, although the correlation between the effect sizes from the two datasets was poor and non-significant. Conclusions This is the first study to comprehensively examine common variants in candidate genes for ASC through meta-analysis. While for majority of the variants, the total sample size was above 500 cases and 500 controls, the total sample size was not large enough to accurately identify common variants that contribute to the aetiology of ASC.